Neutron star properties in a chiral SU(3) model 



M. Hanauske, D. Zschiesche, S. Pai, S. Schramm, H. Stocker, W. Greiner 
Institut fiir Theoretische Physik, J. W. Goethe-Universitat, D-60054 Frankfurt am Main, Germany 



Abstract 

We investigate various properties of neutron star matter within an effective 
chiral SU (3) l x SU (3) r model. The predictions of this model are compared 
with a Walecka-type model. It is demonstrated that the importance of hy- 
peron degrees are strongly depending on the interaction used, even if the 
equation of state near saturation density is nearly the same in both models. 
While the Walecka-type model predicts a strange star core with strangeness 
fraction f s ~ 4/3, the chiral model allows only for f s ~ 1/3 and predicts that 
E , S + and H° will not exist in star, in contrast to the Walecka-type model. 

PACS: 26.60+c, 21.65+f, 24.10Jv 
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Typeset using REVTgX 



The internal constitution and properties of neutron stars chiefly depend on the nature 
of strong interactions. The accepted underlying theory of strong interactions, QCD, is 
however not solvable in the nonperturbative regime. So far numerical solutions of QCD 
on a finite space-time lattice are unable to describe finite nuclei or infinite nuclear matter 
. As an alternative approach several effective models of hadronic interactions have been 
proposed Especially the Walecka model (QHD) and its nonlinear extensions have 

been quite successful and widely used for the description of hadronic matter and finite 
nuclei. These models are relativistic quantum field theories of baryons and mesons, but 
they do not consider essential features of QCD, like approximate SU (3) r x SU(3)l chiral 
symmetry or broken scale invariance. The Nambu-Jona-Lasinio (NJL) model is an 
effective theory which has these features of QCD implemented but it lacks confinement and 
thereby fails to describe finite nuclei and nuclear matter. The chiral SU(3) models, for 
example the linear SU(3)-<7 model [7| have been quite successful in modeling meson-meson 
interactions. The K — N scattering data has been well reproduced using the chiral effective 
SU(3) Lagrangian |||9|]. However, all these chiral models lack the feature of including the 
nucleon-nucleon interaction on the same chiral SU(3) basis and therefore do not provide a 
consistent extrapolation to moderate and high densities relevant to the interior of a neutron 
star. 

This has lead us to construct a QCD-motivated chiral SU(3)l x SU(3)r model as an 
effective theory of strong interactions, which implements the main features of QCD. The 
model has been found to describe reasonably well, the hadronic masses of the various SU (3) 
multiplets, finite nuclei, hypernuclei and excited nuclear matter 0,0] . The basic assump- 



tions in the present chiral model are: (i) The Lagrangian is constructed with respect to the 
nonlinear realization of chiral SU(3)l x SU(3)r symmetry; (ii) The masses of the heavy 
baryons and mesons are generated by spontaneous symmetry breaking; (iii) The masses of 
the pseudoscalar mesons are generated by explicit symmetry breaking, since they are the 
Goldstone bosons of the model; (iv) A QCD-motivated field \ enters, which describes the 
gluon condensate (dilaton field) [12j ; (v) Baryons and mesons are grouped according to their 
quark structure. 

In this letter we investigate the composition and structure of neutron star matter with 
hyperons in this chirally invariant model. The total Lagrangian of the chiral SU(3)l x 
SU(3)r model for neutron star matter can be written in the mean-field approximation as 



(for details see Ref. JTT 
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CSB = - 



Ae P = E Mil^-m^i ■ (2) 

Here £kin is the kinetic energy term of the baryons and the scalar (a, () and vector (ui, (p, p) 
mesons. The interaction Lagrangian of the different baryons with the various spin-0 and 
spin-1 mesons are £bm and £bv, respectively. The sum over i extends over all the charge 
states of the baryon octet (p, n, A, £~, E°, £ + , S - , 5°). £ vcc generates the masses of the 
spin-1 mesons through the interactions with spin-0 mesons, and Cq gives the meson-meson 
interaction term which induce the spontaneous breaking of chiral symmetry. A salient 
feature of the model, the dilaton field Xi which can be identified with the gluon condensate, 
is included. It accounts for the broken scale invariance of QCD at tree level through the 
logarithmic potential. £sb introduces an explicit symmetry breaking of the U(1)a, the 
SU(3)y, and the chiral symmetry. The last term £ lep represents the free lepton Lagrangian. 
The effective masses of the baryons in the nonlinear realization of chiral symmetry are given 

byig 



\g S o^a OS -l){^C-a) 



m N = m 

2 

m* A = m - -goMos ~ l)(v2C ~ a ) 
2 

= m + -go S ( a os ~ l)(v^C ~ °) 
m* s = m + \g s os {2a os + 1){V2( - a) , (3) 

with m = gQ X {\/2a + Q/y/3, in the usual notation |TT[. The parameters g^, gQ 8 and etos 
are used to fit the vacuum baryon masses to their experimental values. The thermodynamic 
potential of the grand canonical ensemble per unit volume at zero temperature for the 
neutron star matter can be written as 

a/v = -£ v „ - £o - c SB -v,-E^/ KM - Kl - i E ^ / -j==- 

(4) 

In Eq. (4) the vacuum energy V vac has been subtraeted. For a given baryon species i, the 
single partiele energy and chemical potential are respectively, 



E*(k) = ^kf + mf , 

Hi = bi/j, n - qiH e = /j,* + g iu] u + g i4> (f)o + g ip hipo , (5) 

with p* = E* (k = fcpj; h and are the baryon number and charge of the ith species. The 
energy density and pressure follows from the Gibbs-Duhem relation, e = Q/V + J2k=n f^kPk 
and P = —VI /V. At a given baryon density ps, the field equations as obtained by extremizing 
fl/V are solved self-consistently in conjunetion with the charge neutrality and /3-equilibrium 
conditions. 
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The parameters of the chirally invariant potential, ko and &2, are used to ensure an 
extremum in the vacuum, while k% is constrained by the rj' mass, and k\ is varied to give 
m a = 500 MeV. The vacuum expectation value of the fields a and ( are constrained by 
the pion and kaon decay constants, a = f n and Co — —{^Jk — fn)/v2- The parameters 
Qj\fij and xo are used to fit the binding energy of nuclear matter B/A = e/ p b — rn N = 
— 16 MeV at the saturation density p = 0.15 fm -3 . In the present calculation we have 
employed the parameter set Cl of Ref. [TT|. The predicted values of effective nucleon mass, 



incompressibility, and symmetry energy at the saturation density are m* N /m^ = 0.61, K = 
276 MeV, and a sym = 40.4 MeV. The remaining couplings to the strange baryons are then 
determined by the additive quark model constraints: 

2 r, V 9'B,<t> ^2 /„\ 

9\lu = 9y,uj = 2g~uj = -Qnlu = 2g 08 ; = 9z<t> = = —9nu ■ (6) 

Figure 1 shows the energy per baryon as a function of baryonic density p b for varying 
neutron-proton asymmetries, 5 = (p n — p p )/pb, calculated in the chiral model. The curve 
5 = describes infinite symmetric nuclear matter with a minimum at po- With increasing 
asymmetry, (5 > 0) the binding energy decreases and the saturation density is shifted to 
lower values. The binding in nuclear matter for small values of S stems from the isospin 
symmetric nuclear forces. At asymmetries 5 > 0.84 (i. e. a neutron to proton ratio > 11), 
the system starts to become unbound even at the low density regimes. The stiffest equation 
of state (EOS) is obtained for pure neutron matter with 5=1. Due to the /3-equilibrium 
conditions, the EOS for neutron star matter (labeled NS) composed of neutrons, protons and 
electrons (npe) is softer as compared to pure neutron matter. The gravitational attraction 
provides the necessary binding of neutron stars. The present results from the chiral model 
corroborate those obtained in Walecka-like models [13,14] and in the relativistic Brueckner- 
Hartree-Fock calculations [ 15 ]. 



Let us now discuss the inclusion of hyperons. With the choice of the parameter set dis- 
cussed above, the chiral model is found to produce unrealistically large hyperon potential 
depths in nuclear matter in comparison to the experimental values of Uj( s ~ —28 MeV for 
the A and S particle. Parameter sets that reproduce reasonable values of Uj^ s are, how- 



ever, found to yield unsatisfactory nuclear properties ||TTJ. Fortunately, explicit symmetry 
breaking can be introduced in the nonlinear realization without affecting, e. g., the partially 
conserved axial-vector currents relations. This allows for the inclusion of additional terms 



for the hyperon-scalar meson coupling [11 



£ hyp = m 3 Tr + ^ [V>, S]) Tr (X - X ) , (7) 

where X represents a scalar and Sg = — [v3(A 8 )^ — <5j]/3 with A's are the usual Gell-Mann 
matrices. In the mean field approximation this leads to the following additional mass term 



m,: = m, 



+ am 3 \V2(5 - S ) + (C - Co)l , (8) 



where m* is given by Eq. (3). With a = n s , where n s is the number of strange quarks in 
the baryon, and, with the parameter m 3 adjusted to = — 28 MeV, the other hyperon 
potentials obtained are = +3.2 MeV, and U§ ~ +30 MeV. Since the potential of S in 
ground state nuclear matter is still not satisfactory, we have used a = 1 as an alternative 
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parametrization. For = —28 MeV, one obtains now = +3.2 MeV, and Uff = —42 
MeV. We are well aware that our choice of the parametrization of the hyperon potentials 
is not unique; the examination of different ways of generating the experimentally favored 
values for and US is in progress and will be reported in a different context. Hereafter 
we have employed a = 1 in our calculations. 

Let us compare the results obtained in the chiral model with that of a Walecka-type 
model. The latter model employed here has a cubic and quartic self-interaction for the a 
field, U(cr) = g2°" 3 /3 + g^a^/A, introduced to get correct nuclear matter compressibility ||, 
and a quartic self-interaction for the vector field u, C v = g^ui^ui^) 2 . This modification 
leads to a reasonable reproduction of the Dirac-Brueckner calculation We have used 
the parameter set TM1 of Ref. fl6|] which gives a saturation density and binding energy 
of p = 0.145 fm" 3 and B/A = -16.3 MeV. It is to be noted that in the TM1 set, the 
values of m* N /m^ = 0.63, K = 281 MeV, and a sym = 36.9 MeV which primarily influence 
the bulk properties of neutron star matter are nearly identical to those in the chiral model. 
The hyperon couplings a-Y are obtained from a depth of Uy = —28 MeV, while the uo-Y 
couplings are obtained from SU(6) symmetry relations (6). Two additional strange mesons, 
( and 0, are introduced which couple only to the hyperons. The (-Y couplings are fixed by 
the condition W w 2U£ w -40 MeV |7gg. 

Figure 2 shows the particle fractions versus the baryonic density in /3-equilibrated matter 
in the chiral model (top panel) and the Walecka-type model, TM1 (bottom panel). With 
increasing density, it is energetically favorable for nucleons at the top of the Fermi sea to 
convert into other baryons. Note that the sequence of appearance of the hyperon species is 
the same in both models. The first strange particle to appear is the S - , since its somewhat 
higher mass (as compared to the A) is compensated by the electro-chemical potential in the 
chemical equilibrium condition (5). Because of its negative charge, charge neutrality can 
be achieved more economically which causes a drop in the lepton fraction. More massive 
and positively charged particles than these appear at higher densities. These are in fact 
generic features found in neutron star calculations with hyperons |p!3| , |14| , |19| , p0| . Because of 
the equilibrium condition (5), the threshold densities of the different hyperon species are, 
however, strongly dependent on the magnitude of the scalar and vector fields at a given 
density and their interactions with the baryons. The nucleon-nucleon interaction which 
determines the variation of neutron chemical potential with density also determines the 
threshold. The relatively smaller attractive scalar fields and thereby larger effective mass 
of all the baryons cause a significant shift of the density at which the hyperons (£°, S + , S°) 
appear in the chiral model calculation. This enhances the saturation values of both the 
neutron and proton fractions (at a level of ~ 30%) in the chiral model. Consequently, the 
cores of massive stars in the TM1 model (and also in other Walecka-type models fT8"| ) are 
in general predicted to be very (> 50%) hyperon-rich. The chiral model, on the other hand, 
predicts neutron stars with considerable smaller hyperon fractions in the core. 

The composition of the star is crucial for neutrino and antineutrino emission which can be 
responsible for the rapid cooling of neutron stars via the URCA process. It was demonstrated 
|2~T|1 that for a npe system rapid cooling by nucleon direct URCA process is allowed when 
the momentum conservation condition fcp + kp e > kp n , corresponding to a proton fraction 
Y p > 0.11, is satisfied. The magnitude of Y p at a given p b, in turn, is determined mainly 
by the symmetry energy, which is nearly identical in the two models. Therefore, in both 
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models this condition is satisfied at densities ps ~ 2.2p thus, rapid cooling by direct URCA 
process can occur. 

The most profound implication of the constitution of neutron star matter on its bulk 
properties are manifested in the equation of state (EOS). The pressure P versus the energy 
density e is displayed in Fig. 3 for neutron star matter in the chiral (thick lines) and TM1 
(thin lines) model. For npe stars (solid lines), although the incompressibility K and the 
effective nucleon mass m* N at the normal nuclear matter density are similar in the two 
models, the EOS at large densities is found to be considerably softer for the chiral model as 
compared to the EOS for TM1. In the latter model, the m* N rapidly decreases with density, 
consequently the EOS passes quickly to one that is dominated by the repulsive vector mesons 
(u and p) leading to a stiffer EOS. This has a strong bearing on the mass and radius of such 
stars discussed below. The structure of static neutron stars can be determined by solving 



the Tolman-Oppenheimer-Volkoff equations p2 |. We use the results of Baym, Pethick and 



Sutherland [Q to describe the crust consisting of leptons and nuclei at the low-density 
{p B < 0.001 fm~ 3 ) EOS. For the mid-density regime (0.001 < p# < 0.08 fm 3 ) the results 
of Negele and Vautherin [E4J are employed. Above this density, the EOS for the relativistic 
models have been adopted. The masses M of the nonrotating neutron star sequence is shown 
in Fig. 4 as a function of central energy density e c in the two models. The corresponding 
mass-radius relationship is presented in Fig. 5. It is observed that the stiffer EOS for a 
npe star in the TM1 model can support a larger maximum mass of M max = 2.16M with 
a corresponding radius of i?M max = H-98 km at a central baryonic density of p c = 6.08po- 
The corresponding values obtained in the chiral model are M max = 1.84M , i?M max = 11.10 
km, and p c = 6.98po- The large difference in the M max values obtained in the two models, 
with nearly identical K and m* N values at po, clearly demonstrates that measurements 
of maximum masses of neutron (npe) stars at high densities cannot be used to constrain 
the incompressibility and the effective nucleon mass around the nuclear matter densities. 
Constraints on K and m* N from radius measurements of massive stars will be even more 
uncertain, since about 40% of these stars' radius originates from the low density EOS. In 
fact, no precise radius measurements currently exist. 

When the hyperon degrees of freedom are included, the EOS (represented by dash- 
dotted lines in Fig. 3) for both the models are appreciably softer as compared to npe 
stars. This is caused by the opening of the hyperon degrees of freedom which relieves some 
of the Fermi pressure of the nucleons. Also the decrease of the pressure exerted by the 
leptons (they are replaced by negatively charged hyperons to maintain charge neutrality) 
contributes to softening of the EOS. Since the threshold density for the appearance of the 
hyperons, especially (S°, S + , H°) (see Fig. 2), are smaller in TM1 model, these stars contain 
more baryon species. This leads to an enhanced softening as compared to that in the 
chiral model. In fact, both models with hyperons predict quite similar values of pressure at 
moderate and high densities; the structures observed in the EOS (Fig. 3) correspond to the 
population of different hyperon species. These should reflect both in the masses and radii 
of the stars. The maximum masses and corresponding radii of stars with hyperons in the 
two models are almost identical with M max = 1.52M , i?M max = 11.64 km, and p c = 5.92p 
in the chiral model, while in the TM1 model, M max = 1.55M , i?M roa3t = 12.14 km, and 
p c = 5.97po- The magnitude of the central densities indicates that the hyperons (E°, S + , H°) 
are entirely precluded in stars for the chiral model whereas all hyperon species appear with 
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comparable abundances in the TM1 model for the maximum-mass star. The strangeness 
fraction, f s = J2i\Si\Pi/ Pb, are vastly different at the center of the maximum-mass stars: 
0.33 vs. 0.75 in the chiral and in the TM1 model, respectively. Thus, models with similar 
nuclear matter incompressibilities and effective nucleon masses, leading to similar maximum 
star masses and corresponding similar radii can however have a widely different baryonic 
constitution! 

For progressively smaller central densities, the masses of the stars with hyperons are 
larger in the TM1 model (see Fig. 4), although the pressure at a given density is smaller 
than in the chiral model. This is because the masses of stars are determined by the overall 
EOS, and the TM1 model possess a distinctively stiffer EOS at moderate densities (see Fig. 
3). The radii of these small mass hyperon-rich stars are quite distinct in these two models 
(see Fig. 5), because the radius of the stars depend most sensitively on the low density 
behavior of the EOS. In both models, stars of mass 1.44M Q (corresponding to the lower 
limit imposed by the larger mass of the binary pulsar PSR 1913+16 also contain many 
hyperon species with sizeable concentration. 

In conclusion, we have investigated the composition and structure of neutron star matter 
in a novel chiral SU (3) model, and compared its predictions with that of a Walecka-type 
model based on different underlying assumptions. The two models with nearly identical val- 
ues of effective nucleon mass, incompressibility, and symmetry energy at the normal nuclear 
density yield widely different maximum neutron star masses and radii. When the hyperon 
degrees of freedom are included, the maximum masses and the corresponding radii in the two 
models are found to be rather similar. However, softness of the nucleonic contribution in the 
chiral model precludes the hyperons S°, S + , S° leading to much smaller hyperon abundances 
in these stars. 
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FIGURES 




FIG. 1. The binding energy per nucleon B/ 'A versus the baryonic density Pb/po for different 
values of neutron-proton asymmetries 5 = {p n — p p )/pb in the chiral model. The normal nuclear 
matter density is po = 0.15 fm -3 . The curve labeled NS describes a neutron star matter consisting 
of nucleons and electrons. 
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Pb/Po 



FIG. 2. The composition of neutron star matter with hyperons in the chiral (top panel) and in 
the TM1 (bottom panel) model. The normal nuclear matter density is po = 0.15 and 0.145 fm -3 
in the chiral and TM1 model. 
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FIG. 3. The equation of state, i. e. the pressure P, versus the energy density e. The results are 
for npe stars (solid lines) , and for stars with further inclusion of hyperons and muons (dash-dotted 
lines). The calculations are in the chiral model (thick lines) and in the TM1 model (thin lines). 
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FIG. 4. The mass as a function of central energy density £ c /eo for npe stars and for stars with 
further inclusion of hyperons and muons in the chiral and TM1 models. The chiral model results are 
represented by thick lines, and the TM1 results are given by the thin lines. The circles correspond 
to the respective maximum masses. The central energy density at normal nuclear matter value is 
£o = 142.9 and 137.0 MeV/fm 3 in the chiral and in the TM1 model, respectively. 
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FIG. 5. The mass as a function of radius for neutron stars consisting of npe and for stars with 
further inclusion of hyperons and muons in the chiral and TM1 model. The symbols have the same 
meaning as in Fig. 4. 
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